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Ultra-high frequency acoustic resonances (^2 GHz) trapped within the glass core (^1 /im diam- 
eter) of a photonic crystal fiber are selectively excited through electrostriction using laser pulses of 
duration 100 ps and energy 500 pJ. Using precisely timed sequences of such driving pulses, we achieve 
coherent control of the acoustic resonances by constructive or destructive interference, demonstrat- 
ing both enhancement and suppression of the vibrations. A sequence of 27 resonantly-timed pulses 
provides a 100-fold increase in the amplitude of the vibrational mode. The results are explained 
and interpreted using a semi-analytical theory, and supported by precise numerical simulations of 
the complex light-matter interaction. 



INTRODUCTION 

The interactions between optical fields and the vibra- 
tional modes of matter have long been a subject of study. 
Over the last two decades, it has become possible to con- 
trol the strength of these vibrations by shaping the tem- 
poral (amplitude and phase) profile of an optical pulse. 
Demonstrations of such "coherent control" have been re- 
ported for example for excitons in quantum wells [1] , THz 
lattice vibrations in crystals [2], photoisomerization of 
retinal molecules [3] multiphoton ionization [4] and very 
recently in an optical microcavity [5]. Along with such 
controlling techniques, the possibility of precisely fabri- 
cating micron-scale structures had raised renewed inter- 
est in the field due to the higher mechanical frequencies 
(> GHz) [6] and stronger susceptibility to optical forces 
of such tiny structures [7], allowing nonlinear effects to 
occur at relatively low optical powers [8]. The purpose 
of this paper is to report coherent control of acoustic 
resonances (ARs) in the small solid glass core of a non- 
linear photonic crystal fiber (PCF) [9]. In these fibers, 
the cladding region surrounding the core is formed from 
an array of parallel air channels, see Fig. l(a,d). The 
large impedance-mismatch between glass and air creates 
strong confinement for acoustic waves incident upon the 
core-cladding interface, and in some cases a phononic 
band gap forms in the periodic cladding. As a result, 
the core (diameter ~1 fim) acts essentially as a cavity 
for transversely-propagating acoustic waves, giving rise 
to trapped ARs with fundamental frequencies in the few- 
GHz range [10]. 

A short optical pulse propagating along an optical fiber 
excites, via electrostriction, acoustic waves [11]. In stan- 
dard (all-solid glass) fibers a large number of ARs are si- 
multaneously excited, giving rise to a complicated acous- 
tic response after propagating several km long fibers [11]. 
In ultra-small-core PCFs the situation is dramatically 
different, however, because tight confinement of acoustic 



(and optical) waves favors excitation of only a few ultra- 
high- frequency ARs localized in the core. The temporal 
acoustic response, which can be detected after only a few 
dozen meters, is then the superposition of a small number 
of damped sinusoidal vibrations with lifetimes inversely 
related to the glass viscosity and acoustic leakage rate. 

Here we report experiments showing that, if a second 
pulse is launched after the first (and well within the life- 
time of the acoustic resonance) , a given AR can be coher- 
ently reinforced or canceled, depending on the time-delay 
between the pulses. Furthermore, using a long sequence 
of pulses equally spaced in time, it is possible to strongly 
excite the particular AR whose frequency matches the 
repetition rate of the pulses 

EXPERIMENT 

In our experiments we used two different PCFs with 
forward Brillouin scattering spectra that had been care- 
fully characterized [10]. The first (PCF#1, 60 m long) 
had a strong AR at 1.96 GHz, with a lifetime of r 11 
ns (mechanical quality factor Q — 135) and the second 
(PCF#2, 10 m long) had one at 1.53 GHz, with a hfe- 
time of r ~ 18 ns (Q=173). Each fiber also displayed 
weaker ARs at higher frequencies: 2.39 GHz for PCF^l 
and 1.83 GHz for PCF#2. 

The experimental set-up used to excite and detect the 
ARs consisted of two external cavity diode lasers, a pump 
laser oscillating at 1530 nm and a CW probe laser at 1550 
nm. An external modulator was used to create trains of 
100 ps (FWHM) pulses with widely controllable inter- 
pulse spacings. The pulses were optically amplified, pro- 
viding pulse energies of roughly 500 pJ (coupled to the 
PCF). After independently controlling the polarization 
state of the two lasers, we co-launched them into the 
PCF using microscope objectives. The polarization axis 
of the probe beam was aligned at 45 ° to the PCF's bire- 
fringent axes, while the pump was polarized along its 
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FIG. 1: Coherent control of acoustic phonons. The cladding pitch and core diameter are 1.67 & 1.27 /im for PCF#1 (tow 
row), and 1.92 & 1.75 fim for PCF#2 (bottom row). The white bars are 1 fim long, (b) Response for a single pulse excitation 
of PCF#1 (blue curve), excitation by two pulses spaced by 260 ps (red curve) and excitation by two pulses spaced by 512 ps 
(green curve), (e) The same as b) but for PCF7^2 and pulse separations of 329 ps for the red curve and 658 ps for the green 
curve. The cartoons in (b) and (e) illustrate the corresponding pump pulse sequence schemes. (c)(f): The Fourier transform 
power spectra corresponding to the time waveforms shown in (b) and (e), corresponding curves share the same color; the arrows 
indicate the multiplication factor used in each curve. 



fast axis (the differential group delay between the pump 
and probe laser light was negligible). At the output, the 
pump laser light was removed using a spectral filter. A 
polarization-controller, followed by a polarizer, acted as 
our analyzer, converting the phonon-induced polarization 
modulation into amplitude modulation, which was then 
detected using a PIN photodiode and a 2.5 GHz real- 
time oscilloscope. The pump pulse duration (100 ps) was 
roughly 5 times shorter than the AR periods (510 ps and 
658 ps) of the two PCFs, so that the excitation can be 
viewed as impulsive. Each pump pulse creates mechani- 
cal strain in the core through electrostriction [12], and by 
adjusting the delay between pulses one can induce either 
constructive or destructive interference in the AR. 

In Fig. l(b,e) we show oscilloscope traces of the inter- 
ference between two sequentially excited ARs for PCF#1 
and PCF#2. By adjusting the delay between the two 
driving pulses we were able to enhance (green curves), 
or suppress (red curves) the amplitude of the main AR 
with respect to single pulse excitation (blue curves). The 
oscillations observed in the red curves in Fig. l(b,e) are 
caused by asynchronous ARs that are not suppressed by 
the second pump pulse, as well as residual oscillations of 
the main AR. The inital spike in all the curves is due to 
the fast electronic Kerr nonlinearity (cross-phase mod- 
ulation) in the fibers. Figures l(c,f) show fast Fourier 
transforms (FFTs) of the time-domain responses in Fig. 



l(b,e), after removing these Kerr spikes. In Fig. Ic, it is 
clear from the red curve that the 1.96 GHz AR has been 
suppressed by a factor of 30 relative to the response for 
a single pulse (blue curve), so that the neighboring AR 
at 2.39 GHz dominates the residual oscillation. PCF#2 
behaves in a very similar manner (Fig. If). These re- 
sults show that one can carry out selective acoustic spec- 
troscopy in such structures by preferentially exciting or 
suppressing targeted ARs. In Supplementary Video 1 we 
show how the different ARs are excited when the inter- 
pulse delay is scanned from 10 to 790 ps. 

To test this technique further, we generated time- 
sequences of pump pulses spaced by the AR period: 510 
ps for PCF#1 and 658 ps for PCF#2. In Fig. 2(a,b) 
we compare the temporal responses of these fibers when 
a single pulse (red curves) and a sequence of 27 pulses 
(black curves) are used. To illustrate the selectivity of 
the process, we show in Fig. 2c the FFT power spectrum 
of PCF^^l under these two excitation conditions. There 
is 100-fold enhancement in the acoustic power of the 1.96 
GHz AR. 

To explore how the AR amplitudes grow with pulse 
count, we plot on Fig. 2d their FFT spectral power as 
a function of the number of pump pulses for PCF#1 
(red circles) and PCF#2 (black triangles). As the pulse 
count increases beyond 20 there is a saturation in both 
cases, even though the time- windows corresponding to 
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FIG. 2: Multi-pulse excitation of ARs. (a) Response of PCF^^l when excited by a single pulse (red curve) and a train of 27 pulses 
(black curve), (b) Same as (a) but for PCF #2. (c) Fourier transform power spectra obtained from the response of PCF#1 to 
a single pulse (red dashed curve) and a sequence of 27 pulses (black solid curve); the arrow indicates the multiplication factor 
used in the red curve, (d) Spectral power amplitude for the target AR as the number of pulses increases. The black and red 
curves represent PCF^^l and PCF7^2 respectively. The dots represent the experimental data, the dashed line the expected 
amplification, and the solid line the expected amplification with saturation parameter S — 0.06. 



these 20 pulses are 9.7 and 12.5 ns for PCF#1 and #2. 
One would expect the saturation to start when the num- 
ber of pulses approaches Nsat ~ t/T {T is the inter- 
pulse period); this yields Ngat — 22 and Ngat = 27 
for PCF#1 and #2. The expected amplification in the 
acoustic power spectrum for a sequence of N pulses is 
(l -|- a; -I- ... -I- a:^~^) , with x = e~^l'^ . This formula is 
plotted in Fig. 2d for both PCF#1 (dashed black line) 
and PCF#2 (dashed red line). In practice however, the 
output pulse energy from the saturated erbium-doped 
fiber amplifier decreases with the number fc of pulses in 
the sequence (approximately as (1 -f kS)~^ , where 5 is a 
saturation parameter) and, therefore, the acoustic power 

increases as {J2k=o i+fcs )^- This formula is also plot- 
ted in Fig. 2d with the best fit value S = 0.06 for both 
PCF#1 (solid black line) and PCF#2 (soHd red hue), in 
good agreement with the experimental data. 



THEORY 

To model the observed coupling to the optical fields, 
one must first solve the elastic wave equation in the pres- 
ence of both viscosity and electrostriction. To accomplish 
this we decomposed the fiber motion u into a superposi- 
tion of all the acoustic modes u = J2m ^™ I"") ^^P '^^mt, 
where is the modal angular frequency. Within the lin- 
ear elastic approximation, one arrives at damped-driven 
harmonic oscillator equations for each Am, which can 
be solved analytically for impulse excitation E (r, i) = 
EqS {t/tQ) \ ^n), where is the pulse duration, Eq is the 
field amphtude, and is the normahzed (('4'„|'5'„) — 
1) electric field vector corresponding to the guided opti- 
cal mode (HEii-like). The driving term is obtained by 
projecting the electrostrictive force on to each acoustic 
mode, i?„j = (u„j|Fs4) / (u,„|u„j). The electrostrictive 
force density is calculated from the electric field vector, 
the elasto-optic constants (pi2 and pm), and the refrac- 



tive index n according to Fst{Tc,t) ~ ^n^[pi2V|E|^ + 
2p44V • (E • E''')] [12]. Because of the radial symmetry of 
the optical mode, when pump light is coupled to one po- 
larization state (say ic), the force term with pi2 will point 
radially (V |E| c>c r) and drive only radial motion. The 
term with 7)44, however, will point along the polarization 
direction (V • (E • EI') cx ic) and drive a motion where 
the core expands in one direction while it contracts in 
the other. 

The change in the refractive index is obtained by 
adding the dielectric tensor perturbation Aem, induced 
by each mode via the elasto-optic effect and weighted by 
the excitation coefficient [12, 13]. The refractive in- 
dex modulation measured by the probe laser coupled to 
the mode and scattered to the mode \^s) is given 

by 

An{t)^ ^ VB„.(*,|Ae,„|*,)X„(t), 

peaCTiefflE^l ^ 

(1) 

where J7p is the pump pulse energy, £9 is the vacuum 
permittivity, «e// is the optical mode effective index, 
c is the velocity of light in vacuum, and p is the sil- 
ica density. The time-dependent function is given by 
^ra(i) — w,7j^ sin(wmt) exp(— rmi)^(t), where d {€) is the 
Heaviside step-function and LOm = \/^1n ~ T^„, where 
Vra is the modal damping parameter. In our experi- 
ment we are interested in inter-polarization scattering 
and, therefore, we choose and I'S's) as orthogonally 
polarized states. 

Focusing our attention on PCF^^l, Fig. 3a shows the 
numerically calculated excitation coefficient for both lin- 
ear (sohd blue circles) and circular (open red circles) 
polarization states of the pump pulses. The T/?2m-like 
modes are not excited with circular polarization while 
the i?om-like ones are polarization insensitive. Therefore 
Fig. 3a is useful to identify the character of the ARs: 
overlapping open red and blue solid circles identify Rom- 
like modes while purely blue circles identify Ti?2m-like 
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FIG. 3: Numerically calculated electrostriction excitation and elasto-optic refractive index modulation for PCF#1. (a) Elec- 
trostrictive excitation for pump light polarized linearly (blue curve) and circularly (red curve); the arrows indicate the most 
optically active modes, shown in parts (c) and (d). (b) Real (red dashed curve) and imaginary (black solid curve) parts of the 
refractive index perturbation spectrum; the inset shows the temporal response, c) l^exx + Ae^y for the i?oi-like mode at 2.318 
GHz. d) Ae^jy for the ri?2i-like mode at 1.956 GHz. The black solid lines represent the undeformed fiber geometry used in 
the calculation. The mesh deformation in the transverse plane was scaled for clarity. 



modes. Due to a dense clustering of modes near = 0, 
we have plotted only those with Bm > 5%. The strongest 
ARs in this figure (indicated by arrows) have frequencies 
of 1.956 GHz (Ti?2i-like) and 2.318 GHz (i?oi-like). The 
profiles of displacement and dielectric tensor perturba- 
tion are shown in Fig. 3(c,d) where the radially symmet- 
ric (3c) and torsional-radial (3d) motion can be observed 
(animated versions are available in the Supplementary 
Videos 2 and 3). The mode amplitudes after a excitation 
by single 500 pJ pulse are rather small, 4 fm typically, 
due to the high frequency of the motions. 

In Fig. 3b we show the calculated frequency de- 
pendence of the refractive index modulation (assuming 
Tm = 9.1 X 10'' s~^). The red dashed and black solid 
curves correspond to the real and imaginary parts. The 
corresponding time waveform is shown in the inset. The 
index modulation spectrum shows a strong peak only at 
1.956 GHz (Ti?2i-like mode), with much weaker mod- 
ulations near 2.32 GHz, in contrast to the excitation 
coefliicients (Fig. 3a), where several modes contribute 
significantly. This is because most of the excited modes 
(large \B,n\) do not efficiently perturbs the dielectric ten- 
sor (small |Aem|). On the other hand, modes that the- 
oretically have almost no contribution to An might be 
observed in the experiments. This happens, for example, 
for the i?oi-like mode at 2.318 GHz. In a perfectly sym- 
metric PCF (such as the one used in the simulations), 
the i?om-like family behaves much like in a standard op- 
tical fiber in the sense that it does not induce any signifi- 
cant birefringence and therefore does not couple together 
orthogonally polarized optical fields. This results in a 
very small anisotropic index modulation for such modes. 
However, in real PCFs, small geometric distortions will 
slightly alter the exact symmetry of the i?om modes, al- 
lowing them to induce a degree of birefringence which 
can be detected in our experimental set-up. Indeed the 
i?oi-like mode identified in Fig. 3a oscillates at 2.318 
GHz, which is close to the experimentally- measured AR 



at 2.39 GHz. 

In conclusions, ARs trapped in the small core of highly 
nonlinear PCFs, oscillating at GHz frequencies, can be 
coherently controlled by sequences of optical driving 
pulses with durations shorter than a half-cycle of acous- 
tic oscillation. Such multi-pulse excitation can be used 
to preferentially amplify or attenuate selected resonances. 
Control of ultra-high frequency resonances may have ap- 
plications in mode-locking of fiber lasers at very high 
repetition rates, in the realization of optically- pumped 
acoustic "sasers" , or in the design of more efficient fiber- 
based acousto-optic devices. 

We would like to acknowledge N. Joly for fabricating 
the fibers, and the financial support from the Brazilian 
agencies CAPES, CNPq, and FAPESP. 



[1] A. P. Heberle, J. J. Baumberg, and K. Kohler, Phys. Rev. 

Lett. 75, 2598 (1995). 
[2] T. Feurer, J. C. Vaughan, and K. A. Nelson, Science 299, 

374 (2003). 

[3] V. I. Prokhorenko et al, Science 313, 1257 (2006). 
[4] A. V. Sokolov et al., Phys. Rev. Lett. 8703 (2001). 
[5] N. D. Lanzillotti-Kimura et al., Phys. Rev. Lett. 99, 
217405 (2007). 

[6] T. Carmon and K. J. Vahala, Phys. Rev. Lett. 98, 
123901 (2007). 

[7] M. Eichenfield et al, Nat. Photonics 1, 416 (2007). 

[8] T. Carmon et al, Phys. Rev. Lett. 94, 223902 (2005). 

[9] P. St. J Russell, J. Lightwave Technol. 24, 4729 (2006). 
[10] P. Dainese et al.. Opt. Express 14, 4141 (2006). 
[11] E. M. Dianov et al. Opt. Lett. 15, 314 (1990). 
[12] A. S. Biryukov, M. E. Sukharev, and E. M. Dianov, 

Quantum Electronics 32, 765 (2002). 
[13] D. Marcuse, Bell Syst. Tech. J. pp. 985-995 (1975). 



